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FOREWORD 


This report has been prepared to expedite early dissemination of the 
information generated under the contract. The data and conclusions must be 
considered preliminary and subject to change as further progress is made on 
this program. This is a progress report covering the work done during the 
first 12 months of the contract, and it is not a final report. 
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1.0 


INTRODUCTION 


Modern jet engine design imposes extremely high loadings and temperatures 
on hot section components. Fuel costs dictate that minimum weight components 
be used wherever possible. In order to satisfy these two criteria, designers 
are turning toward improved materials and innovative designs. Along with 
these approaches, they also must have more accurate, more economical, and more 
comprehensive analytical methods. 

Numerous analytical methods are available that can, in principle, handle 
any problem that might arise. However, the time and expense required to pro- 
duce acceptable solutions is often excessive. This program addresses this 
problem by setting out a plan to create specialized software packages which 
will provide the necessary answers in an efficient, user-oriented, stream- 
lined fashion. Separate component-specific models will be created for burner 
liners, turbine blades, and turbine vanes using fundamental data from many 
technical areas. The methods developed will be simple to execute, but they 
will not be simple in concept. The problem is extremely complex and only by a 
thorough understanding of the details can the important technical approaches 
be extracted. The packaging of these interdisciplinary approaches into a 
total system must conform to 'he modular requirements for useful computer pro- 
grams. 


The overall objective of this program is to develop and verify a series 
of interdisciplinary modeling and analysis techniques that have been special- 
ized to address three specific hot section components. These techniques will 
incorporate data as well as theoretical methods from many diverse areas 
including cycle and performance analysis, heat transfer analysis, linear and 
nonlinear stress analysis, and mission analysis. Building on the proven tech- 
niques already available in these fields, the new methods developed through 
this contract will be integrated to provide an accurate, efficient, and uni- 
fied approach to analyzing combustor burner liners, hollow air-cooled turbine 
blades, and air-cooled turbine vanes. For these components, the methods 
developed will predict temperature, deformation, stress, and strain histories 
throughout a complete flight mission. 

This program, to a great extent, draws on prior experience. This base of 
experience is invaluable for understanding the highly complex interactions 
among all the different technical disciplines as well as for estimating the 
importance of different engine parameters. In particular, there are four spe- 
cific areas in which experience is especially beneficial. 

First, with the recent increases in fuel costs, greater emphasis has been 
placed on more accurate solutions for stresses and strains in order to under- 
stand and improve the durability and life of hot section components. Conven- 
tional linear elastic analyses are no longer sufficient; instead, they now 
provide the boundary values for more refined creep and plasticity calculations 
These nonlinear analyses are now performed routinely as part of the design 
process at General Electric. This extensive experience with these plasticity 
and creep methods contributes directly to developing component specific models 



Second, advances in 3-D modeling capability are being achieved by the 
concepts developed under the NASA-supported ESMOSS pros-am. ESMOSS concepts 
provide the basis to develop an efficient modeling system for geometric and 
discretized models of engine components. 

Third, the NASA-funded Burner Liner Therma l/St t uc tural Load Modeling Pro- 
gram contributes strong support to this program. The specific area addressed, 
transfer of data from a 3-D heat transfer analysis model to a 3-D stress analy- 
sis model, will provide the background and framework for the data interpolation 
required for all thermomechanical models in this contract. 

Fourth, over the past 10 years General Electric has developed internally 
a family of computer programs: LASTS, OPSEV, and HOTSAM. These programs ail 

have the common thread of using selected points from cycle data, heat transfer, 
and stress analyses and a decompos itior./synthesis approach to produce accurate 
values of temperature, stress, and strain throughout a mission. These programs 
are totally consistent with the overall objectives of this program, and repre- 
sent a proven technology base upon which the component specific models are 
being developed. Significant advances being made are the inclusion of non- 
linear effects and the introduction of improved modeling and data transfer 
techniques . 

The program is organized into nine tasks which can logically be separated 
into two broadly parallel activities (Figure 1). On the right of Figure 1 we 
have the Component Specific Thermomechanical Load Mission Modeling path. 

Along this path a Decompos i t ion/ Synthe s i s approach is being taken. In broad 
terms, methods are being developed to generate approximate numerical models 
for the engine cycle and the aerodynamic and heat transfer analyses needed to 
provide the input conditions for hot parts stress and life analysis. 

The left path, Component Specific Structural Modeling provides the tools 
to develop and analyze finite element nonlinear stress analysis models of com- 
bustor liners and turbine blades and vanes. These two paths are shown in more 
detail in Figures 2 and 3. 

Software Development, Task IV, consists of planning and writing the com- 
puter programs for both paths, with the necessary interconnections, using a 
structured, top down approach. 

In the Thermomechanical Load Mission Modeling portion of the program 
(Figure 2), we are developing in Task III a Thermodynamic Engine Model which 
generates the engine internal flow variables for any point on the operating 
mission. The method for doing this is described below. Task V is developing 
techniques to decompose flight missions into characteristic mission segments. 

In Task VII a Thermomechanical Mission Model is being developed. This uses 
the flow variables from the Thermodynamic Model to determine metal temperature 
and pressure distributions for a representative combustor liner and turbine 
blade and vane. 
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Figure 2. Component Specific Thermooechanical Load 
Mission Modeling,. 
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Figure 3. Component Specific Structural Modeling. 
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Individual Casks for the Structural Modeling activity are shown in Figure 
5. The requirements of Software Design, Task II, have been factored into Task 
fl, the evaluation of the structural analysis mechods which were selected for 
•valuation in Task I. Task VIII provides the capability for structurally 
aodeling current state-of-the-art combustor liners and hollow turbine blaaes 
ind vanes, given the defining dimensional parameters. These parameters will 
>e chosen to facilitate parametric studies. 

The component specific models are being developed in two steps. In the 
first a geometric model is defined. In the application of the Component 
Specific Modeling Program these data are then transferred to che Thermomeehan- 
.cal Load Mission Model to provide the geometry for determining component 
>ressures and temperatures. Thus, a data transfer link is being developed to 
Io this in Task IV, Software Development. The capability for generacing from 
: he geometric model a discretized, finite element model is also a part of Task 
fill. At this point another link between the two paths is needed to transfer 
:he component temperatures and pressures from the Thermomechanical Load Model 
:c Che finite element model, interpolating the data as needed to define nodal 
: emperatures and pressures. This also is being completed in Task IV. 


o 


2.0 TECHNICAL PROGRESS 


2. 1 TASK I - LITERATURE SURVEY 

The first task of this program was to perform a literature survey of 
available methods, techniques, and solution strategies that can be used to 
geometrically model, display, and structurally analyze burner liners, tur- 
bine blades, and vanes. NTIS, NASA, DTIC, and internal General Electric 
Company documents were searched. The fruits of this survey are listed in 
Appendix A. 


2.2 TASK II - DESIGN OF STRUCTURE ANALYSIS SOFTWARE ARCHITECTURE 


The software architecture was designed us-ing the methodology developed 
on the ESMOSS program. The first step in this process was to perform a 
functional analysis of the problem using Soft Tech's Structured Analysis 
and Design Technique (SADT). This analysis defined the work of the software 
system and provided the foundation for the software development of Task 
IV. This analysis was performed by teams whose members had expertise in 
all of the pertinent areas. Appendix B contains the diagiams showing the 
functional decomposition and the current data dictionary. 

Figure 4 shows the software architecture defined for the structural anal- 
ysis portion of COSMOS. This is the architecture developed for internal pro- 
grams. Subroutine ELKIND determines the element type, the output option, the 
material number, the orthotropic orientation set, the stiffness computation 
code, the area load set, and the line load set. This information is then bit 
packed into one word. Subroutine NODDOF determines the degree of freedom per 
node consistent with the element type. Subroutine FIXITY determines zero dis- 
placement boundary conditions and sets counters to eliminate the proper 
equations from the solution. Subroutine PRETIS determines all other prescribed 
boundary conditions. Subroutine CONNEC reads the element nodal connectivity. 
Subroutine CONSTR establishes any constraint equations. Subroutine XYZCOR 
reads the global coordinates cf each node. Subroutine MIDNOD generates 
midside nodes for those elements requiring them. Subroutine SKWDOF reads any 
skew boundary conditions. Subroutine MTA5LE reads material property data. 
Subroutine ELDATA establishes the element specific data. Subroutine COLHT 
establishes the column heights for the linear solver subroutine. Subroutine 
SBLOCK determines the number of solution blocks required by the linear equa- 
tion solver. Subroutine FILCOR allows for more portability of of code bv 
establishing logic to determine whether there is sufficient core or whether 
files are needed. Subroutine ELSTIF develops the element stiffness. Sub- 
routine ELMASS develops the element mass. Subroutine ELDAMP develops the 
element damping. Subroutine ADDSTF assembles the global stiffness matrix. 
Subroutine MODSTF modifies the global stiffness matrix for boundary conditions. 
Subroutine EIGEN performs an eigenvalue solution and outputs the answers. 
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Subroutine LOADS develops the right-hand side vector. Subroutine MODLOD modi- 
fies the right hand side vector for the prescribed boundary conditions. Sub- 
routine SOLVER solves the set of linear equations. Subroutine NONCON develops 
the nonlinear constitutive equation information. Subroutine NONSTF develops 
the nonlinar stiffness terms for a tangent modulus solution methods. Subrou- 
tine NONRHS develops the nonlinear right hand side terms for a right hand side 
pseudoforce solution method. Subroutine OUTPUT develops the requested 
results and sends them to the proper I/O device. This structure is both 
modular and growable. 


2.3 TASK III - THERMODYNAMIC ENGINE MODEL 


The Thermodynamic Engine Model has been completed. The model has been 
developed as a simple calculations tool which will take as inputs the three 
variables altitude *.h) , Mach number (M) , and power level (PL;or the allowed 
flight map of an engine, as shown in Figure 5. In addition, ambient tempera- 
ture deviations from the standard atmosphere, airframe bleed air requirements, 
and engine deterioration can also be included as- part of the input to the 
thermodynamic model. For each input condition specified Ly h, M, and FL, the 
thermodynamic model will calculate gas weight flow (w) , temperature (t), and 
pressure (p) at selected aerodynamic engine stations as needed to determine 
component thermal loadings. These stations are shown in Figure 6. 

The technique for developing a thermodynamic engine model is shown in 
Figures 7 and 8. The engine to be analyzed must be defined thermodynamically 
by an engine cycle deck (computer program) that can be run to generate the 
internal flow variables at the chosen aerodynamic stations (Figure 7). To 
encompass the complete engine operating map (Figure 5), 148 operating points 
are chosen and w, t, and p are calculated using the cycle deck for the 
selected stations as weli as N^ and N 2 , the fan and core speeds. From this 
station data, an engine performance cycle map is constructed. This is essen- 
tially a set of three-dimensional data arrays that map the station data (w, 
t, p, N^, and N 2 ) onto the engine operating map (Figure 5). Given an arbitary 
operating point defined by h, M, and PL, it is then possible in principle to 
interpolate on the engine performance cycle map to determine station data. In 
practice the station parameters are nonlinear functions of the input param- 
eters, and considerable effort is needed to develop these multidimensional 
interpolations. The computer programs used to generate the engine performance 
cycle map from the engine cycle deck output has been developed as part of Task 
III. The functioning of the thermodynamic engine model is shown in Figure 8. 
Given an engine mission, as shown schematically in Figure 9, it can be defined 
by values of the input variables h, M, and PL at selected times through the 
mission. Using these input variables and the engine performance cycle map, an 
interpolation program developed in this effort will calculate engine station 
parameters throughout the mission (Figure 8). These are then used to define 
station mission profiles of w, t, p, N[ , and N 2 as functions of time at each 
aerodynamic station. These station mission profiles then become the input to 
the thermomechanical engine model. 
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Figure 5. Engine Operating Map. 
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Figure 6. CF6 Aerodynamic Stations. 
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2.3.1 


Detailed Specification and Requirements 


The first step in this task was to develop detailed specifications and 
requirements for the thermodynamic engine model software. This specification, 
which was approved by the NASA Program Manager, is presented in Appendix C. 

2.3.2 Model Design and Development 

Based on the detailed specifications, the thermodynamic engine model 
software was generated and a set of 148 performance cases was obtained to 
load this model. The next task was to establish interparameter interpolation 
functions. To assist in this effort, 103 special cases that would maximize 
interpolation errors were chosen from the cycle deck. 

A master interpararaeter linearity study was executed to evaluate inter- 
polation functions. A computer program (STATPAC) was available that could 
take 30 input performance parameters, perform transformations on the data, 
and generate crossplots of the transformed function. The linearity of these 
crossplots was the criterion of excellence in the selection of interpolation 
functions. One hundred validation cases were run with 30 parameters each, 
giving 3000 individual comparisons. Sixty-three additional performance 
cases were used to perform the interparameter linearity study for t..e Mach 
number and altitude control variables. 

Based on the above program, a set of interparameter interpolation func- 
tions and transformation functions was defined and encoded in the TDE model 
software. The accuracy that can be achieved with these is excellent. As a 
final "trial' run," this model was tested against the CFM56 engine flight con- 
ditions . 

The TDE user's manual is in preparation. 

2.4 SUPPORTING TASKS 

2.4.1 Task IV - Software Development 

This task consists of planning and writing the computer codes for 
paths of this program with the necessary interconnectors. As such, it 
continuous and ongoing effort; the substance is covered under the othe 
headings . 

2.4.2 Task VI - Structural Analysis Methods Evaluation 

The objective of Task VI was to evaluate the basic methods to be 
developing the structural-analysis capabilities of the component-spec^ 
models. The task has been completed, and selected items are being imp 
in Task VIII. 


both 
is a 
r task 


used in 
f ic 

lemented 


12 



Burner Cas Average Discharge Temperature, 


Ch 


Design 



Vane Overall Cooling Effectiveness, 


H = (T, - T , , , 

4 VA n S 


) / ( T - T ) 
X' v 4 2 


Figure 10. 


Turbine Vane Cooling Effectiveness. 
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2.5 TASK V - MISSION MODEL DEVELOPMENT 


The thermomechanical, model is still under development. The form is based 
on types of correlations previously developed within GE . Figure 10 shows a 
representative correlation for a turbine vane. Metal temperatures at various 
points on the vane T are correlated in terms of a vane overall cooling 
effectiveness n, station gas temperature Tj at compressor discharge, and T 4 
at combustor discharge. Using the station mission profiles, it is possible to 
calculate the temperatures at selected locations on each component as func- 
tions at time, given the input parameters h, M, and PL that define the engine 
mission. These provide the boundary conditions for the component stress and 
analyses . 


2.5.1 Combustor Liner Temperature and Pressure Decomposition 
and Synthesis 

Introduction 


Work has been completed on a generalized procedure established to predict 
liner temperatures and pressure drop for a rolled ring combustor. This pro- 
cedure was developed using available data, both measured and and calculated, 
from several sources. The correlating procedure was demonstrated using engine 
test data at sea level conditions. Since no engine data at altitude conditions 
were available, high pressure sector data for a CF6-80A rolled ring liner were 
used to demonstrate that the method could be used at altitude conditions. A 
THTD (transient heat transfer - Version D) • two-dimensional calculation for the 
rolled ring combustor at takeoff conditlrn was lsed to generate the cooling 
effectiveness data for the entire panel (Panel / outer) that was selected for 
analysis. The development of the procedure and the steps leading to the final 
coolir.g effectiveness curve and pressure drop data are described below. 


Liner Temperature Distribution 

A large amount of engine test data is available for the CFb-51) rolled 
ring combustor. For the test of the engine, all panels were instrumented 
with the largest amount of thermocouples on Panels 3 and 7 of the outer liner. 
Using the metal temperature data, the combustor exit temperature T 4 , and the 
compressor discharge temperature T 3 , a cooling effectiveness n c was evaluated 
from. 


^4 ~ ^liner 

T4 ' 


(1) 


The results for Panel 7 outer are shown in Figure 11 for several condi- 
tions from idle to takeoff power. As shown in the figure, the method corre- 
lates the data for several instrumented locations over a wide range of sea 
level conditions and also shows a slignt increase in cooling effectiveness 
for several of the liner postions as the power level (pressure) decreases. 
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This pressure effect is illustrated in Figure 12 where Che ratio of the effec- 
tiveness at a given power to the effectiveness at takeoff power is plotted 
versus pressure. As shown, the ratio for Panel 7 outer increases to about 
1.03 at the idle condition. This ratio was plotted for all panels. It was 
found that the ratio increased to about l.Ott on tne forward panels and was 
about 1.03 for all panels downstream of the aftmost dilution panel. This dif- 
ference is probably due to the difference in the contributions of flame radi- 
ation and convection to the total heat load to the liner. The heat load on 
the forward panels is controlled by radiation, which is a power function of 
the local gas temperature. The gas temperatures are lower and Che gas 
velocities are higher on the aft panels. Thus, the heat load is controlled 
by convection and is a direct function of the local gas temperature. If the 
pressure effect is not used, then Che predicted temperatures at the idle con- 
dition would be only about 13* F higher on the aft panels than if the correc- 
tion is used. The difference on the forward panels would be about 55* F at 
the idle conditions. Therefore, the pressure correction can be neglected and 
not cause serious errors in the stress/life calculations. 

Data at altitude conditions were not available from the engine test; 
however, both altitude and sea level data were available .'rom a high pressure 
sector test of a rolled ring combustor, which has the same type of liner con- 
struction as the engine test combustor. These data, converted to cooling 
effectiveness, are snown in Figure 13 and illustrate that the generalized 
procedure can be used at both altitude and sea level conditions. 

As mentioned earlier, a THTD analysis of Panel 7 was used to generate a 
cooling effectiveness curve for the entire panel length. This curve was gen- 
erated for maximum (hot streak) and average panel temperatures. Generally, hot 
streaks exist around the circumference of the liners and streak for each fuel 
injector. These streaks result from locally high gas temperatures and veloc- 
ities near the liner surface. Data matching calculations have led to analyt- 
ical procedures to simulate hot streaks in the THTD predictions. These pro- 
cedures used factors applied to local gas temperature and velocity to generate 
the hot-streak heat transfer input. These curves are shown in Figure 14 and 
the coordinates of the curves for input to the computer model are given in 
Table I. The coordinate system used in TrlTU analysis and the model node 
layout are shown in Figure 15. 

Temperature Gradient Through Material Thickness 

ure gradient through the material thickne 
veness, compressor discharge temperature 
ature. The temperature gradient through 


( 2 ) 


ure, * F 
ture, * F 


ss 

and 

the 


An expression for the temperat 
can be derived from cooling effecti 
pressure, and combustor exit temper 
material can be calculated from 


Tu - T, 


Q/A 

k 


where 


hot side metal temperat 
cold side mecal tempera 
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Table I. Cooling Effectiveness Coordinates for Panel 7 Outer. 



Coordinates , 
Inch 

Cooling Effectiveness 

Locat Lon 

X 

Y 

Hot Streak 

Average 

Liner 

0.034 

0.03 

0.957 

0.957 


0.094 

0.03 

0.925 

0.925 


0.158 

0.03 

0.880 

0.880 


0.208 

0.03 

0.849 

0.849 


0.254 

0.03 

0.832 

0.832 


0.288 

0.03 

0.819 

0.819 


0.340 

0.03 

0.798 

0.802 


0.404 

0.03 

0.773 

0.788 


0.438 

0.03 

0.706 

0.760 


0.584 

0.04 

0.659 

0.739 


0.654 

0.04 

0.632 

0.722 


0.764 

0 . 04 

0.623 

0.708 


0.854 

0.04 

0.645 

0.739 


0.934 

0.044 

0.677 

0.750 


1. 04 

0.044 

0.713 

0.760 


1.074 

0.04 

0.735 

0.767 


1.114 

0.03 

0.740 

0.774 


1.150 

0.03 

0.744 

0.777 

Forward 

0.144 

0.250 

0.968 

0.968 

Cooling 

0.194 

0.240 

0.952 

0.952 

Ring 

0.234 

0.230 

0.945 

0.945 


0.248 

0.210 

0.933 

0.933 


0.284 

0.170 

0.900 

0.900 


0.300 

0.124 

0.865 

0.865 


0.308 

0.084 

0.838 

0.838 

Aft 

0.980 

0.116 

0.791 

0.820 

Cool ing 

0.984 

0.176 

0.881 

0.890 

Ring 

0.984 

0. 234 

0.940 

0.940 


0.794 

0.290 

0.956 

0.956 


1.024 

0.329 

0.967 

0.967 


1.054 

0.34 

0.970 

0.970 


1.090 

0.354 

0.973 

0.973 


1.124 

0.360 

0.982 

0.982 


1.158 

0.366 

0.976 

0.976 


1.204 

0.370 

0.944 

0.344 













Q/A = heat flux through material, Btu/hr-ft 2 
t = material thickness, ft 

k = metal conductivity, Btu-f t/f t 2 -hr-° F 

The heat flux can be estimated from: 

Q/A = (h c + h r ) (T c - T 3 ) (3) 

and is proportional to (Tuner - 

Q/A - (h c + h r ) T L i ner " t 3 > (4) 

where 

hc = convection heat transfer coefficient Btu/hr-ft 2 -* F 
hr = equivalent heat transfer coefficient for radiation to the casing 
Btu/hr-ft 2 -' F 

Substituting the heat flux expression into the gradient Equation (2j gives: 

t H " T c (n c + h r> fe N 

- T - * 

^Liner *-3 k 


Using the equation for cooling effectiveness, 

” ■‘■Liner 
n c = — = 

t 4 “ T 3 

an equation for (Twiner " T 3 ) can be written as follows: 

(T L iner ~ T 3 ) - (1 - h c ) T 4 + (h c - 1) T 3 

Substituting Equation (7) into the expression gives ( 6 ) 

X H ~ T C a (He n r } T 

(1 - n c ) t 4 + (n c - l) r 3 a k 


(7) 


(8) 


The convection term n c varies with pressure; thus, the gradient through 
the material thickness should be correlated with pressure. 

A THTD analysis was done at several pressure conditions and the calculated 
temperature gradients were plotted versus P 3 for several axial locations. 

The results are shown in Figure 16. The locations are indicated in Figure 15. 
As shown in the figure, the gradient data are correlated with pressure. The 
constants m and b in the equation: 


- T, 


l Liner “ T 3 


Tr - T c 

(1 - n c ) T 4 + U c - 1) T 3 


mP 3 ♦ b 


(V) 


are tabulated in Table II. 
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Pressure, psia 


Figure 16. Material Thickness Temperature Gradient. 
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Table II. Linear Fit Constants for Equation ( 8 ). 


Locat ion 

X, inches 

m 

b 

1 

0.094 

12.3 x IO - 5 

0.100 

2 

0.438 

14.1 x IO - 5 

0.061 

3 

0.654 

9.0 x 10-5 

0.061 

4 

0.854 

10.7 x 10-5 

0.092 

5 

1.114 

28.1 x IQ ' 5 

0.168 


Given the combustor exit temperature T 4 , the compressor discharge pres- 
sure P 3 , and the compressor discharge temperature T 3 , the temperature gradient 
through material thickness can be calculated from Equation (9) using che cool- 
ing effectiveness from Table I and the constants from Table II. 

A generalized procedure has been established to predict liner temperature 
distribution and the ':emperature drop through the material thickness. The 
predictions can be maae using the above constants and equations and the engine 
cycle data; that is, the compressor discharge temperature and pressure and the 
combustor exit temperature. 


General Piessure Drop 

The pressure drop data from the engine test were reviewed and com- 
pared to values predicted by the COBRA (Combustor Analysis) program. The 
static pressure tap locations are shown in Figure 17 and the comparisons of 
predicted and measured pressure drop data are shown in Figures 18, 19, and 20 
as a function of the combustor flow function squared, ( w Comb / PT3 T 3 . The 
subscripts shown on the pressure drop curves refer to the static pressure tup 
locations identified on Figure 17. 

In general, the agreement between the measured and predicted values is 
good except for the values for inner liner. The engine pressure drop to the 
inner passage is larger than the predicted values for both the forward and 
aft locations. Data from a full scale, full annular diffuser test were 
examined and were found to agree with the predicted values. Other data 
from a combustor -test have shown that the forward passage pressures are 
sensitive to the axial location and that more realistic passage pressures are 
obtained if the pressure taps are mounted to the casing wall rather than to 
the combustor cowl or to a film sic'' as was done in the engine test. If the 
inner passage measured data are corrected based on the diffuser test data, 
then the agreement with the COBRA prediction is much better. 

The above comparisons indicate that a COBRA analysis can be used to 
predict the pressure drop across the dome and each of the liners and is 
related to the square of the combustor flow function by: 
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Combustor Flow Function Squared M 'omb 











Symbols Represent Engine Data 



Combustor Flow Function Squared 







AP 

?3 


K 


^Comb 

P T3 


2 


*3 


( 10 ) 


where K is a constant for each of the various combustor components. The 
values of K for the locations shown in Figures 18, 19, and 20 are given 
in Table III. 


Table III. Pressure Drop Constants for Equation (10). 


Locat ion 

7 1 

1 

Dome 

10.0 

X 

10-3 

Outer Liner 




Forward 

9.1 

X 

10-3 

Aft 

13.5 

X 

10-3 

Inner Liner 

9.6 

X 

10-3 

Forward 

15.0 

X 

10-3 

Aft 





Given the combustor flow and the compressor iischarge temperature and pressure, 
the various combustor pressure drops can be calculated. 

The above comparisons were based on data obtained at sea level conditions. 
Additional data were needed at altitude conditions to determine if the sea 
level prediction curves could be used over the entire flight map. Both sea 
level and altitude test data are available from a high pressure test of a 
60° sector of a rolled ring comoustor. The overall pressure drop data are 
shown in Figure 21. The data from both test conditions correlate together as 
shown in the figure. Also shown in Figure 21 is the design pressure 
drop as a function of the sector flow function squared. The COBRA analysis 
can thus be used to predict the combustor pressure drop and can be applied 
over the entire flight map. 


Summary 


Correlations to predict liner temperature, temperature drop through the 
mate .al thickness, and liner pressure drop have been defined. These correla- 
tions use the compressor discharge temperature and pressure as well as the 
combustor exit temperature to calculate the above values over the entire flight 
map. The information needed to develop the computer program is summarized 
below. 
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O Sea Level Conditions 
f~~) Altitude Conditions 





1. Complete axial temperature distributions are calculated from Equa- 
tion (1) using the cooling effectiveness values given in Table I. 

2. Temperature gradients through the material thickness are calculated 
from Equation (9) using the values given in Tables I and II. 

3. Liner and dome pressure drops are calculated from Equation (10) 
using the constants in Table III. 


2.5.2 Turbine Blade and Vane Temperature and Pressure 
Decomposition and Synthesis 

Our objective was to establish a generalized procedure for predicting 
the airfoil metal temperature distribution throughout the engine operating 
power range. Three approaches were initially assessed. They were the follow- 
ing : 

1. One-Dimensional Heat Transfer Balance - A generalized one-d imens ionai 
heat transfer equation was considered for calculating the metal 
temperature distribution across the turbine vane or blade. The 
equation takes the general form 

* f (location, pressure ratio, ( 11 ) 

temperature ratio, gas coolant weight- 
flow ratio, etc) 

2. Predicted Cooling Effectiveness - Use of THTD-pred icted design point 
temperature distributions for vanes or blades to establish the 
generalized cooling effectiveness equation. These design tempera- 
ture predictions were collected from CF 6-6 and CF6-50 design groups. 

A total of 20 local temperatures were selected for each vane or 
blade to establish the local generalized cooling effectivenesses. 

The airfoil metal temperature could then be calculated by the 
equat ion 

T m ■ Tg - n c (T g - T 3 ) (12) 

where n c denotes the generalized cooling effectiveness derived from 
the predicted design temperature distributions, and would be a func- 
tion of Tg, T 3 , and airfoil axial location. 

3. Tested Cooling Effectiveness - Similar approach as described in 
Item 2, but with the cooling effectiveness determined from tempera- 
ture measured in engine tests rather than predicted values. 

Measured airfoil temperature data were collected from the CF 6 
single-shank Stage 1 blade and the ''■tage 1 vane. 


T _ 
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These three approaches in predicting the airfoil metal temperatures dis- 
tribution were assessed by further detailed studies described below. 

Engine test data were collected from two current development engines for 
the Stage 1 vane and from previous engine tests for the Stage 1 blade. Using 
the measured temperature along the airfoil surface, T m , the turbine inlet 
gas temperature, T 41 , and the compressor discharge temperature, T 3 , the cool- 
ing effectiveness p c can be calculated from 

T 4 L ~ T m (13) 

" c ■ *41 - t 3 

for each thermocouple location. 

Figure 22 shows the thermocouple location along the airfoil surface for 
the vane and blade. Figure 23 is a plot of the cooling effectiveness p c 
defined in Equation (l) and derived from the engine cest data for the Stage 1 
blade. Beginning at a lower power level and progressing to design power level 
the test results furnish a series of data points that defines the trends of p c 
throughout the engine operating range. Similar plots for the turbine vanes, 
as derived from engine test results, are given in Figures 24, 25 and 26 for 
the vane 152, 502, and 85% span. These nearly linear p c lines derived from 
test data throughout the power level range were used as one of the three 
approaches to synthesize turbine vane and blade temperatures at selected ioca- 
t ions . 

The method using THTD-predicted airfoil metal temperatures from design 
analysis to establish the generalized cooling effect ivenes » equation for the 
vanes and blade is described and summarized in the following. Predicted THTD 
metal temperature distributions for the Stage 1 vane and blade were obtained 
from two current development engines. Figures 27 and 28 show the surface 
metal temperature distribution for the airfoil pitchline at design conditions 
for the Stage 1 vane and blade respectively. Using the predicted temperature 
along the airfoil surface T m , and the turbine inlet gas temperature T 4 [, and 
the compressor discharge temperature T 3 , the generalized cooling effectiveness 
p c can be evaluated from Equation (13). T^iese generalized p c ' s based on 
predicted metal temperatures are plotc:c on Figures 29 and 30 for the Stage 1 
vane and blade, respectively, compared with the corresponding p c deduced from 
engine test data. As can be seen from these figures, the predicted p (denoted 
by solid symbols) and the p c extended from test data (open symbols) have the 
same relative relationship to each other, but absolute quantities do not 
match. Similar plots are shown in Figures 31 and 32 for the Stage 1 vane 15% 
and 852 span respectively. The disagreement between these two p c ' s is not 
surprising, and the reasons are thought to be as follows: 

• Predicted metal temperatures are based on nominal values of design 
variables such as Tg as , T c 00 1 a n t , manufacturing tolerances, and 
thermal properties 
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Figure 22. Thermocouple Location for Vane and blade. 
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Figure 26. Field Plot of Generalized Turbine 
Vane Cooling Effectiveness at 85% 
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Figure 29. Stai?e 1 Vane Test Results at 50% Span. 
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Turbine Inlet Temperature. 



Figure 30. Stage 1 Blade Te?t Results. 
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• Uncertainties in calculating the gas-side and coolant-side heat 
transfer coefficients, the radiation heat flux, film cooling effec- 
tiveness, etc. 

• Uncertainit ies in thermocouple measurements, flow checking measure- 
ments, etc. 


No effort was spent to investigate these uncertainit ies or the design 
variable tolerance or to refine the test data matches for ..he vane and blade. 
THTD is a well-established analytical too. to predict the airfoil metal tem- 
peratures. Rather, our total effort was used to search the THTD predictions 
for blades and vanes at engine idle or intermediate power setting conditions. 

To establish a general procedure to predict the static gas pressure dis- 
tributions along the airfoil surface, typical design gas pressure distribu- 
tions for the engine vane and blade were collected and plotted in Figures 33 
and 34 respectively. These pressure distributions were assessed to develop the 
thermomechanical load model. 

The final general procedure for predicting turbine vane and blade metal 
temperatures and their gas-side static pressure distributions throughout the 
engine operating range is as follows: local distributions of the generalized 

cooling effect ivenesf for the vane and blade are given in Figures 35 and 36, 
and were obtained from the THTD-predicted metal temperatures reported above. 
These cooling-effectiveness distributions are plotted against the normalized 
airfoil surface length S/L, and are shown in Figures 37 and 38. By specifying 
the engine cycle data T 3 and T^i and the airfoil suction and pressure-side 
surface lengths, the airfoil surface metal temperature distribution can be 
gotten from Figures 37 and 38 for the typical turbine vane and blade used in 
this study. 

In a similar fashion, the normalized airfoil pressure distributions for 
turbine vane and blade are given in Figures 39 and 40. By specifying the 
engine cycle data (T 4 for the vane, for the blade and airfoil geometry), 

gas-side static pressure distributions along the airfoil surfaces can be 
seen in these figures. 

The final method to gain a generalized cooling effectiveness from one- 
dimensional heat balance analysis yields the resulting equation: 





0.2 




(14) 


f 

where the * sign denotes the reference condition such as the THTD design point 
or the test data point, and 
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Figure 34. C. as Static Pressure Distribution Along Blade Surface. 











figure 38. Turbine Blade Normalized Cooling Effectiveness Distribution at Pitch Line 






Figure 39. Turbine Vane Normalized Cas Static Pressure Distribution at Pitch Line 





Figure 40. Turbine Blade Normalized Gas Static Pressure Distribution at Pitch Line 




Wc/Wg 

P T e 
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8 


- j , Tg ■ T 4 for vane 
* ^ for blade 

■ vane or blade cooling flow, X Wg 
gas-side total pressure, P ?4 or Pi 4 i 
gas-side temperature, T 4 or T 4 L 


K 


T g - 


T3 


1 


Equation (14) can be used to predict the generalized cooling effectiveness n c 
at any power setting from the referenced n c * and the stage variables Tg, T 3 , 
and the cooling flow (Wc/Wg). The referenced n e * can either be gotten from 
THTD design prediction or test data or simply by the following equation through 
one-dimensional heat balance analysis, 

n c * * nf* ♦ w c G pc 

h g A g 


film cooling effectiveness 
coolant specific heat 
gas-side heat transfer coefficient 
gas-side heat transfer area 
Equation (11) can be restated as follows for local conditions: 


where 


if 


Cpc 


n g 

Ag 


T g - 


t 3 


- 1 


(15) 


n 


c,n 



(16) 


where 

Tg * Gas temperature relative to the surface, * R 

T m * Metal temperature at a specific location, * R 

T c ,i * Inlet coolant temperature, * R 
n c,n 3 Cooling effectiveness at location n 
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The value of n c , n has been shown to vary with power load in an approximately 
linear manner based on thermocouple temperature measurements. Some 
variances have also been noted between based on thermocouple temperature 

measurements and calculated temperatures from detailed 20 or 3D thermal analy- 
sis models. Such differences may result from any of the following causes - 
taken alone or in combination: 

• Thermocouple variability and deterioration 

• Calculated versus actual distribution of gas side heat transfer 
coefficients 

• Unknown true temperature of the gas or coolant in the vicinity of 
the thermocouple 

• Thermocouple is located in a high temperature gradient region 

• 2D model versus 3D actual airfoil geometry 

Use of the cooling effectiveness allows the calculation of a local 

metal temperature, ^m,n, when given the corresponding values of T g , T c> i > 
which are derived from cycle data and n c>n determined for the given airfo’l 
design and operating conditions. The metal temperature is 

Tm , n * ^g ” he , n ( ^g “ ^c , i ) (17) 

It would be convenient if expressions could be derived for predicting nc,n 
without running a detailed thermal analysis for each operating condition 
for stress analysis, as there can never be enough metal thermocouples to 
define n c ,n f° r locations needed by mechanical designers. But it is 
generally possible to install a few thermocouples on airfoils, sufficient to 
refine or verify procedures for calculating airfoil temperature distributions. 
The refinements result from matching the calculated temperatures at test point 
operating conditions which correspond to those for the measured temperature 
data. 


The use of thermal analvsis models, tuned up by temperature matching 
analyses, offers the wealth of detailed temperature distributions needed for 
stress analysis of cooled airfoils. However, it is inefficient and costly to 
run such detailed thermal analysis models for the many power and environmental 
conditions encountered in a flight mission. A cost-effective approach is to 
identify the regions in an airfoil where critical stresses occur. Tempera- 
tures at selected points in these stress regions can be converted to corre- 
sponding values of n c ,n at other power or mission points of interest fcom those 
derived from a detailed thermal analysis model at a specified reference set of 
operating conditions. Such conditions are available at the design point and 
may be available at other conditions. In principle, ne,n can be obtained from 
Equation (18). 
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Based on a 

Detailed Thermo Analysis 
Mode 


( 18 ) 


n c ,n 


n c ,n 


( n c,n,Ref) 


c , n , Re f 

where the subscript Ref stands for a specific set of reference conditions. 


Previous work furnished some preliminary and approximate equations for 
(n c n /n c n Ref^- Further work was done in search of more appropriate approxi- 
mate expressions. These expressions were based on the assumption of a series 
of one-dimensional heat exchangers, dominated by convective heat fluxes, to 
represent the surface of an airfoil. Some variances may result from using an 
approximate expression to relate n c to operating parameters instead of run- 
ning a complete thermal analysis mode.* The relationship of n cn to n c n,Ref as 
determined from two detailed THTD runs of an engine high pressure turbine 
Stage 1 blade and vane at sea level takeoff and altitude cruise conditions was 
studied. 


The approximate equation relating n c>n to n Cini R e f can be tested by sub- 
stituting turbine operat ig parameters for cruise and takeoff conditions and 
comparing the results with n Cjn and n Cjnj R e f based on the detailed THTD ther- 
mal models. The resulting first-order equation is 



P Tg 

P Tc 


Average absolute total pressure of gas stream 
Average absolute total pressure of coolant 


( FF) g 
(FF) C 


t Tc 
T T g 
A e,g 
A f ,c 


w /RT T , 

— a — - — “■ turbine flow function 

A f,g P Tg 
w c^Tc 

= — turbine flow function 

A f,c p Tc 

Average absolute total temperature of coolant, # R 
Average absolute total temperature of gas, * R 
Gas flow area 

Coolant supply, effective flow area for the same airfoil, 
turbine, and engine 
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Af g/Af )C ■ (Af ,g/Af ,c/Ref and f° r normal operating ranges ( FF)g/( FF)g t Re f 
and (FF5 C /(FF) C )Re f are essentially constant. For convenience, Tjg is some- 
times replaced with related temperatures like T 4.1 and T-jg and Tt c by T 3 , 

where 

T 4 .I ■ Average absolute total temperature ahead of the turbine 
rotor, * R 

TXB * Gas temperature relative to the surface, * R 

T 3 * Compressor discharge absolute average temperature, * R 

The use of these engine cycle parameters in Equations (19) results in a 
single first-order expression for the nozzle vane and blade. 


1 - 

1 “ ^c.Ref 


_ \ 0 . 04 

l! ) 

T 3 ,Re f / 


/ t 4 . 1 , Re f \ 

\ T 4.1 ) 


0.04 


For the CF6-50C high pressure turbine: 


T 4 . l 

Takeoff (Reference Conditions) 1555 3078 

Cruise Conditions 4 1317 2446 


from which, 


~ n c , cruise 
1 " He, T/0 


1.00255 


( 21 ) 


THTD 2D calculations of pitch-line vane and blade temperature distributions 
are listed in Table IV. The sea level takeoff data were used as the reference 
conditions. Similar fltitude cruise data were used for the test cases. 

Figure 41 compares (n C)tl ) cruise (n Cifl )x/Q based on Tg and surface tem- 
perature. For the turbine Stage 1 vanes, the trend is similar for the 14 
sampled temperatures identified in Figure 41. The same data are replotted in 
Figure 42 for comparison with the results of Equation (20). It can be seen 
that with the exclusion of three points in the leading edge region (Numbers 1, 

2, and 14) the data group has about an average value of 0.92 as compared to a 
value of unity from Equation (21). Also note that n c n cruise “ n c n T/0 
when (1 - n c , n )/( 1 - n C) n iRe f) - 1. 
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Figure 42. Prediction of Cooling Effectiveness Relationship at Cruise 
to Takeoff From the First Order Equation Versus Detailed 
THTD Results. 
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Table IV. High Pressure Turbine Stage 1 Vane and Blade Temperature 
Distribution. 

• Sea Level T 3 - 1095* F T4.1 - 2618* F 

• Cruise T 3 - 857* F T4.1 - 1986* F 


Vane 

Blade 


Takeoff 

Cruise 


Takeoff 

Cruise 

Location 

^m, n 

^c ,n 

^■ra,n 

n c ,n 

Locat ion 

, n n c , n 

^m, n 

n o,n 

n 

• 

F 

• 

F 

n 

* F 

O 

F 

l 

• 1850 

0.647 

1496 

0.605 

1 

1897 0.473 

1347 

0.566 

2 

1785 

0.677 

1364 

0.687 

2 

1857 0.500 

1341 

0.511 

3 

1884 

0.631 

1387 

0.672 

3 

1863 0.496 

1336 

0.576 

4 

1529 

0.797 

1160 

0.813 

4 

1913 0.463 

1354 

0.560 

5 

1683 

0.725 

1258 

0.752 

5 

1917 0.460 

1364 

0.551 

6 

1759 

0.689 

1309 

0.721 

6 

1923 0.456 

1366 

0.549 

7 

1822 

0.66 

1364 

0.687 

7 

1884 0.482 

1370 

0.546 

8 

1862 

0.641 

1402 

0.663 

8 

1975 0.422 

1382 

0.535 

9 

. 1865 

0.639 

1412 

0.657 

9 

1912 0.464 

1375 

0.541 

10 

1745 

0.696 

1318 

0.715 

10 

1797 0.539 

1326 

0.585 

ll 

1718 

0.708 

1241 

0.763 

11 

1805 0.534 

1340 

0.572 

12 

1633 

0.748 

1227 

0.771 

12 

1821 0.523 

1342 

0.570 

13 

1866 

0.639 



13 

1809 0.53L 

1334 

0.578 

14 

1778 

0.681 

1498 

0.604 

• 14 

1811 0.530 

1328 . 

0.583 

Bulk 

1601 

0.763 

1237 

0.765 

15 

1738 0.578 

1312 

U.597 






16 

1846 0.507 


0.581 






Bulk 

1759 0.564 

1335 

0.577 


T 8 , # * 

» 3233 

2475 


T- ' F - 3233 

6 » 

2475 

T 4 .I. ’ F 

■= 2618 

1986 

T 4 

* F - 2615 

1986 






T TB , ’ F - 2560 

1940 


*3. ' E 

- 1095 

857 


T 3 , * F - 1095 

857 


The foregoing results were anticipated because the vanes are hollow shells 
with impingement inserts. Therefore, the heat source and sink were essen- 
tially constant over the vane surfaces. Obviously, the simple flat plate heat 
exchanger model was better applied to the pressure and suction surfaces of the 
airfoil than to the "cylinder like" leading edge region. 

Figure 43 compares (n C|n ) to (nc, n )T/0 based on T4 . 1 and surface tem- 
peratures for the H? turbine Stage l blade. The data for the 16 points around 
the airfoil lie in a band along a straight line. However, the cooling effec- 
tiveness shown for surface temperatures at cruise power is consistently higher 
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Comparison of Cruise and Sea Level Takeoff Cooling 
Effectiveness . 





Chan at takeoff, with deviations up to 27% near the trailing edge. A somewhat 
better agreement would have been attained if local mean wall temperatures had 
been used. On closer examination, some of the 16 locations sampled are 
affected by 2D conduction in contrast to the ID flat plate heat exchanger 
assumed in the derivation of Equation (4) and (5). The variances can also be 
due to local variation in coolant temperature as the air traverses through 
serpentine cooling passages, and to the variation in film cooling effective- 
ness. Hence, Figure 43 illustrates a fundamental problem in trying to derive 
a simple, but general, expression for predicting airfoil temperature distribu- 
tions when the cooling system design is complex and several modes of heat 
transfer are combined. The data of Figure 43 constitute a correlation, but 
this is based on a detailed design analysis of the particular airfoil with 
cruise boundary conditiors in the thermal model rather than a simple projec- 
tion from the reference takeoff case. This problem is likely to be a common 
occurrence with sophisticated cooled airfoil designs. It is therefore planned 
to make a similar comparison of cruise and takeoff cooling effectiveness of 
another high pressure turbine blade design. 

2.5.3 Stress-Strain Decomposition ana Synthesis 

The decomposition and synthesis of stresses, strains, and deformations is 
technically the most, challenging portion of this program. It requires innova- 
tive methods to produce usable results for burner liners, turbine blades, and 
vanes. Thus, our goal under this task has been to compile a library of possi- 
ble decomposition and synthesis techniques and to assess their validity. 

Among the technioues being considered are the following: 

• Assume that the structure remains totally elastic at all stress 
levels and do the decomposition and synthesis based on an elastic 
"pseudostress . " 

• Assume that the structure is deformation-controlled strain range 
invariance). The first level of decomposition and synthesis would 
be based on deformations (total strains). A second level of syn- 
thesis could then introduce the effects of plasticity and creep by 
using the material response characteristics to partition the total 
strain into elastic, plastic, and creep components. 

• Assume that the structure is load-controlled (stress invariance). 
Decomposition and a first level of synthesis would be based on load 
terms reflecting the centrifugal loadings and the temperature and 
pressure distributions. A second level of synthesis could then 
introduce the effects of plasticity and creep by using the material 
response characteristics to determine the elastic, plastic, and 
creep strains that would be caused by the total load. 

• Use simplified nonlinear finite element modeling to decompose and 
synthesize the stresses, strains, and deformations in terms of the 
set of analyzed mission components. These simplified models could 
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be either one 2D or 3D element or a nonlinear substructure. These 
models could use boundary conditions from the detailed analysis or 
they could be run as an intimate part of the detailed analysis. 

• Apply the method ot superposition for the decomposition and synthe- 
sis of stresses, strains, and deformations. This method would be 
investigated based on the following hierarchy of calculated parame- 
ters : 

- deformations 

- strains 

- stresses 

We will determine to what degree these parameters can be decomposed 
and synthesized by superposing the results from individual loading 
functions (temperature, pressure, rpm) . 

• Use linear and nonlinear interpolation of the results of a detailed 
analysis for decomposing and synthesizing stresses, strains, and 
deformations. The interpolating parameters would be second-level 
predicted temperatures, pressures, and rpm's. 

• Form look-up tables of deformations, stresses, and strains as func- 
tions of temperatures, pressures, and rpm's. These tables would 
then be used to decompose and synthesize the mission cycles. 

• Finally, generate from test data an empirical model relating 
stresses, strains, and deformations to temperatures, pressures, and 
rpm's. With 'this model, mission cycles could be decomposed and new 
ones synthesized. 

Our first step was to survey existing techniques tor decomposition ana 
synthesis. The best documented or these was used in combustor design. Com- 
bustor life is limited by (1) creep of the outer liner, (2) oxidation, (3) high 
cycle fatigue, and (4) low cycle thermal fatigue. Thermal fatigue is the most 
limiting on the CF6 family of engines. Life projections are made with a tool 
called CO-LIFE analysis. In this, thermal strains are computed as 

c^x * aAT 

Metal temperatures are computed using a fiim effectiveness 

^me t a 1 ~ ^c 
^gas ” Tc 


The mechanical pressure strains are computed as a linear function of the pres- 
sure differential 

e A p 3 CAP 
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The CO-LIFE analysis Cool is used Co predicc insealLed life. As such, ic 
musC evaluace derace condicions and make facCory Co field life projeccions. 

To do Chis ic utilizes acCual engine paramecers and actual test results to 
develop the necessary constants. The internal logic includes the following 
calculat ions : 

Ac - (ea ♦ e^p) 

Eat = K 1 ( c HOT “ c COLD^ 

3 k 2 (t HOT ■ TcOLd) 
e AP ■ k 3 p 3 

t HOT * t 3 + n c (T 41 “ T 3> 

Ac =• KN f “ + K 5 N 4 0 

Laboratory 100-second hold time smooth bars test data are used to define the 
material life capability. 

This calculation method was subsequently upgraded to the HOTSAM program 
by G. Weber, one of the COSMOS team. The following changes were made in the 
CO-LIFE logic. 

• • Use (chot “ c COLD) instead of (T^ot " T COLD^ 

• Assume K varies with metal stiffness so that 


6 AT = K 1 (e HOT " e COLD ) i e C0LD- 

Our goal in COSMOS is to improve these techniques for decomposing and 
synthesizing the temperatures and pressures and the stress-strain response. 

Two well-documented hot section problems are used as the tools in the 
initial efforts at decomposition and synthesis. The first of these is the 
uniaxial model from NASA CR-165268, "Turbine Blade Tip Durability Analysis." 
This model simulates the strain-temperature-time conditions occurring at a tur- 
bine blade squealer-tip cracking location. As a multiaxial tool we are using 
the problem reported in NASA CR-2271, "Multiaxiil Cycle Thermoplasticity Anal- 
ysis with Besseling's Subvolume Method." The problem analyzed in this paper 
is that of a shingled combustor segment with a hot spot. Table V gives the 
pertinent information for the critical location of the squealer tip problem. 
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The simplifying approach of assuming strain range invariance for this 
problem is obviously feasible and involves only a 2.8% error. However, the 
other individual parameters vary too much to make this approach acceptable by 
itself. 


Table V. Results of Turbine Blade Tip Analyses, Inelastic . 


• 

Elastic 

Cyc le 
1 

Cyc le 
2 

Cyc le 
3 

Maximum Total Strain, % 

0.025 

-0.05 

-0.062 

-0.0829 

Minimum Total Strain, % 

-0.2925 

-0.3582 

-0.371 

0.3918 

Total Strain Range, % 

0.3175 

0.3082 

0.3090 

0.3089 

Mean Stress, ksi 

-23.9 

-1.9 

4.9 

' 11.7 


Figure 44 shows the temperature-time cycle, and Figure 45 gives the total 
strain-tine cycle, which are the imposed boundary conditions for the problem. 
Using the material properties as given in NASA CR-165268 and the classical 
constitutive equation embedded in CYANIDE, Figures 46', 47-, and 48 show the 
results of the computer prediction. Figure 46 is the stress-time response, 
Figure 47 is the plastic strain-time response, and Figure 48 is the creep- 
time response. This will be used as the baseline against which decomposition/ 
synthesis techniques will be measured. 

Figure 49 shows the stress versus time prediction if the problem is con- 
sidered totally elastic. Figure 50 shows the stress versus time response and 
Figure 51 shows the creep strain versus time response when plasticity is 
ignored . 

Table VI is a comparison of the pertinent data for the three types of 
problem simulations. As can he seen, an elastic analysis (E) gives no mean- 
ingful data, whereas the approximation ignoring plasticity (EC) is useful. 

The mean stresses and stress ranges approximate reality. 
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Table VI. Comparison Data for Three Simulated Problems. 


Temperature, * F 



1 

2 

3 

4 

5 

6 

n 

8 

E 

EC 

EPC 

650 

650 

650 

1900 

1900 

1900 

1800 

1800 

1800 





n 



; 1 

in. ) 

-2550 

-2550 

-2550 

-2925 

-2925 

-2925 

1 

I 

-600 

-600 

-600 

0.0 

0.0 

0.0 

• R - St 

S 

EC 

EPC 

ress (ks 

0.0 

0.0 

0.0 

i) 

-54.4 

-50.4 

-30.1 

-53.5 

-49.0 

-26.8 

-49.0 

-36.2 

-23.4 

-49.5 

-21.8 

-16.4 

7.0 

50.6 

55.4 

-14.8 

23.9 

28.1 

0.0 

43.6 

48.4 

• R - ?1 

E 

EC 

EPC 

as tic Sc 

0.0 

0.0 

0.0 

— 
rain ( ir 

0.0 

0.0 

-1222 

.. /in. ) 

0.0 

0.0 

-1222 

0.0 

0.0 

-1262 

0.0 

0.0 

-1262 

0.0 

0.0 

-1129 

0.0 

0.0 

-1129 

0.0 

0.0 

-1129 

• R - Cr 
E 

EC 

EPC 

eep Sere 

0.0 

0.0 

- 

lin ( in. 

0.0 

-211 

lJ!L_ 

in. ) 

0.0 

-217 

-52 



0.0 

-762 

-255 

i 

0.0 

-1568 

-611 

0.0 

-1568 

-611 

0.0 

-1568 

-611 



0.0 

-1568 

-611 




Legend 


E = Elastic only 
EC * Elastic and Creep 
EPC * Elastic, Plastic, and Creep 


Next, the turbine blade tip durability model was exercised with a simpli- 
fied thermal cycle. Reference is made to Figures 44 and 45 for the original 
complex cycle being investigated. To assess the effect and utility, the tem- 
perature cycle was revised as shown in Figure 52. Table VII shows the 2D 
CYANIDE results for this revised cycle. These values compared to those in 
Table VI show this approximation to be very good. 
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Table VII. Simplified Temperature Cycle Results. 





Temperature 

, * F 




1 

2 

3 

4 

5 

6 

7 

8 

E 

650 

1900 


1900 

1900 

600 

600 

600 

EC 

650 

1900 

1900 

1900 

1900 

600 

600 

600 

EP 

650 

1900 

1900 

1900 

1900 

600 

600 

600 

EPC 

650 

190U 

1900 

1900 

1900 

600 

600 

600 

• R - Tot 

al Str) 

tin 







E 

0.0 

-2856 

-2550 

-2925 

-2800 

256 

-600 

0.0 

EC 

0.0 

-2850 

-2550 

-2925 

-2800 

250 

-600 

0.0 

EP 

0.0 

-2850 

-2550 

-2925 

-2800 

250 

-600 

0.0 

EPC 

0.0 

-2850 

-2550 

-2925 

-2800 

250 

-60U 

0.0 

1 

• R - Stress 



• 





E 

0 

-54.4 

-48.7 

-55.9 

-53.5 

7.0 

-16.7 

0.0 

EC 

0 

-50.4 

-44. f 

-39.1 

-23.1 

51.3 

27.6 

44.3 

EP 

0 

-31.1 

-25.4 

-31.6 

-29.2 

42.4 

18.7 . 

35.5 

EPC 

0 

-30.1 

-24.4 

-27.2 

-18.6 

55.0 

31.4 

48.1 

• R - P 1 a 

is tic S 

train 



• 




E 

0 

0 

0 

0 

0 

o 

0 

0 

EC 

0 

• 0 

0 

0 

0 

0 

o 

•o 

EP 

0 

-1222 

-1222 

-1273 

-1273 

-1273 

-1273 

-1273 

EPC 

0 

-1222 

-1222 

-1238 

-.1238 

-1136 

41136 

-1136 

1 

• R - Creep Str 

ain 







E 

o 

0 

0 

0 

0 

0 

0 

0 

EC 

0 

-211 

-215 

-877 

-1591 

-1591 

-1591 

-1591 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

-51 

-52 

-263 

-590 

-590 

-590 

-590 


In the multiaxial work, the model under analysis is a single shingle com- 
bustor segment analyzed as a flat plate in a condition of plane stress. The 
shingle segment is shown in Figure 53 and modeled as illustrated in Figure 54. 
The combustor shingle model was chosen because of its multiaxial nature and 
the complex thermal cycle. The thermal condition of the combustor shingle at 
peak temperature is shown in Figure 55. Figure 56 defines the thermal cycle 
at the center of the hot spot. 
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Figure 55. Temperature Distribution on Shingle 
at Peak Condition. 
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Figure 56. Thermal Cycle at Center of Hot Spot. 
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Figure 57. Effective Stress Versus Effective 
Strain at Center of Hot Spot for 
Base Case (1) . 
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The combustor model was run under four separate analysis conditions to 
enable the effects of creep and plasticity to be separated. The analysis con 
ditions were as follows: 


• Elastic/Plastic/Creep 

• Elastic/Plastic 

• Elastic/Creep 

• Elastic only 


(HSR80) 

(HSR80-B) 

(HSR80-C) 

(HSR80-D) 


Figures 57 through 71 give pertinent information concerning the state of 
stress for the element at the center of the hot spot. 


2.6 TASK VIII ~ COMPONENT SPECIFIC MODEL DEVELOPMENT 
2.6.1 Geometric Modeling 

• 

The geometric modeling of these three specific components will dovetail 
with and make use of the modeling capability being developed in the ESMOSS 
contract. At the end of the ESMOSS contract, the personnel involved will be 
joining this program, bringing with them the new modeling capability. Mean- 
hile, all of the necessary hooks to other parts of this computer system are 
now being generated. 


2.6.2 Remeshing and Mesh Refinement 

This area and the next, Self-Adaptive Solution strategies, touch on each 
other synergistical ly . What is sought in this program is the best' combination 
of both. This involves two major areas of investigation: the method to be 

used to refine, upgrade, and rearrange the mesh, and the criteria to be used to 
activate this process. 

There are a number of ways to refine a mesh to get a better answer: 

(1) one way is to progressively subdivide a coarse mesh, always retaining all 
previous meshes within the finer mesh; (2) a second family of techniques 
totally realigns the mesh based on some criteria such as strain 
energy density; (3) a third method is to leave the mesh unchanged but upgrade 
the order of the elements. 

The first method, progressive subdivision, has certain theoretical and 
computational advantages. If the finite element interpolating functions used 
meet the requirements for completeness and continuity, convergence is mathe- 
matically guaranteed when we refine the mesh by progressive subdivision. The 
computational process of remeshing by progressive subdivision is straightfor- 
ward; however, it guarantees a larger problem to solve. 
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Figure 58. R-Direct^on Stress-Strain Cycle at 
Center of Hot Spot for Base Case (1) 
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Figure 59. Z-Direction Stress-Strain Cycle at 
Center of Hot Spot for Base Case(l) 
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Figure 60. Effective Stress Versus Effective 
Strain at Center of riot Spot -for 
Case 2. 
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Figure 61. R-Direction Stress-Strain Cycle at 
Center of Hot Spot for Case 2. 
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Figure 64. Effective Stress Versus Effective 
Strain at Center of Hot Spot for 
Case 3. 
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Figure 65. R-Direction S tress- S train Cycle at 
Center of Hot Spot for Case 3. 
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Figure 66. Z-Direction Stress-Strain Cycle at 
Center of Hot Spot for Case 3. 



Figure 67. R-Direction Stress Versus Z-Direction 

Stress at Center of Hot Spot for Case 3. 
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Figure 68. Effective Stress Versus Effective 
Strain at Center of Hot Spot for 
Case 4. 
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Figure. 70. Z-Direction Stress-Strain Cycle at 
Center of Hot Spot for Case 4. 
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Figure 71. 

R-Direction Stress Versus Z-Direction 
Stress at Center of Hot Spot for Case 4 







For a solution of the finite element system of equations: 

[K] (4) • (F) 

suppose there is a numerical solution for the displacement, {6*}. Then the 
equilibrium or residual force vector is generated: 

{R} - {F} - [K] {«*} 


A perfect solution would result in this vector containing all zeros. 

Given the finite numerical accuracy of the computer, this is impossible. 
Therefore, a measure of the numerical "goodness" of the solution is to be 
found in how much this vector deviates from zero. Decisions on whether to 
re-solve or redefine the problem can be based on the total and local devia- 
tions from zero. If a few local degrees of freedom are out of equilibrium, 
this might suggest a local remeshing. If the total equilibrium is deficient, 
this will require remeshing and/or re-solving with greater numerical accuracy. 

The decision tree for this is as follows 


1 . 


If IRi < C R 
and all R^ < Cr^ 


the solution is good 


2 . 


If IRi > Cr 

and [Number of nodes with Cr^ < Ri 
then re-solve 


Cr- 


iu 


> C 


s 


3. If IRi > Cr 

and [Number of nodes with Ri > Cr. ] < C s 

then remesh and re-solve 


4. If Ri < Cr 

but some Ri > Cr. 

then remesh and re-solve 


where : 

Ri ■ i c ^ residual-free vector 

Cr ■ Maximum allowable sum of Ri 

Cr^ * Lower bound for Ri for possible remeshing 

Cr^ * Maximum allowable upper bound for Ri 

Once an acceptable d isj'l acement solution has been reached, proceed to the 
element level. If, at the elastic level, stresses and strains are linearly 
connected, only one of these two needs to be evaluated. Strain will be 
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checked. The total strain at each calculation point in an element is made up 
of an elastic strain and a thermal strain: 


One aspect of this program is the establishment of acceptable strain 
gradients for different element types. Between adjacent strain calculation 
points in one element, and probelly over the entire element, a strain gradient 
would not he chosen that could encompass an elastic-plastic-elastic or a plas- 
tic-elastic-plastic variation. Therefore, 

if 


«* - c * > I 2‘yield 1 • 
i J 7 

remesh this element. 

Additionally, there will be a change in sign. Therefore, 
if 



remesh this element. 

Once the nonlinear solution has been entered, the element level checks 
become more complex and more important. The total strain is now made up of 
the elastic strain, thermal strain, plastic strain, and creep strain: 


i i i i i 

Now stress and strain are no longer linearly connected; stress is a 
function of elastic strain only. Once again, between any two adjacent calcu- 
lation points within one element, an elastic strain gradient greater than the 
allowable material elastic gradient is not desirable. Thus, 


if 


C . £ . ^ I Id | , 

1 J 


remesh this element. The limit on the thermal strain would still be retained. 
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If 


c _ 

— — < 0 f remesb this element. 

c 

j 

The next check is on the computed plastic and creep strain. No sign 
changes in either of these are allowed. In addition, a maximum gradient is 
set . 

if 



j 

or 



j 

or 

I c° - c° | > C 
i j P 

or 



> C , 
C 


remesh this element. 


Next, proceed to the interelement level check. These are of the same 
nature as the above, but now involve adjacent calculation points in adjacent 
elements. 


2.6.3 Se lf-Adaptive Solution Strategies 

In the development of basic self-adaptive solution strategies, we are 
using the work of Edward T. Wilson of the University of California at Berkeley, 
and Joseph Padovan and Surapong Tovichakehaikul of the University of Akron. 
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WiLson's efforts are directed toward an overall solution strategy, while 
Padovan's and Tovichakchaikul 1 s work is on load incrementing and time-stepping 
for geometrical and material nonlinear solutions. 

Wilson's philosophy on internal program organization for SAP-80 computer 
programs is applicable to the Component-Specific Modeling Program, with some 
extensions. He suggests that the basic internal organization of a computer 
program for structural analysis depends strongly on the method used to form 
and solve linear equations, with the frontal and profile (or active column) 
methods most often used. Both have the exact same economy so that the choice 
must be based on other factors. 

In the frontal method, element stiffnesses and solutions of equations are 
formulated in a joint sequence manner. Therefore, all element stiffness sub- 
routines, the equation solver, and the front of the stiffness matrix must be 
in core storage (or rolled in and out) during the reduction of the stiffness 
matrix. For the profile approach, the formation of all element stiffnesses 
for a particular type of element can be accomplished by a single call to one 
program segment. The formation of the total stiffness is a separate program 
segment in which the element stiffnesses are read in sequence from secondary 
storage and the total stiffness matrix is formed in active column blocks. In 
this case, the actual solution phase is another separate program link. Evalu- 
ation of substructure stiffnesses, calculations of mode shapes and frequencies, 
and evaluations of reactions and member forces are all sepacate links. This 
clear uncoupling of different phases of the program gives the profile approach 
a clear advantage in modularity and adaptive solution techniques. Also, the 
profile approach has no significant disadvantages when compared to the frontal 
me thod . 

Padovan and his coworkers at the University of Akron have been developing 
"Self-Adaptive Incremental Newton-Raphson Algorithms" for nonlinear problems. 
They use a three-level approach. In the first Level, incremental Newton- 
Raphson operators are used to "tunnel" into the problem solution space. Tne 
second level involves the constant monitoring of the different stages of solu- 
tion through various quality/ccnvergence/nonl inearity tests. The third 
level works with the results of the second level. The violation of any of the 
quality/convergences/nonlinearity tests triggers various scenarios for modify- 
ing the incremental Newton-Raphson strategy. The self-adaptive modifications 
triggered by the third level fall into one of three categories: global stiff- 

ness reformation; preferential, local reformation; or load increment adjust- 
ment. Recently, they have developed constrained, self-adaptive solution pro- 
cedures for structures subject to high temperature elastic/plastic/creep 
effects. In this, they used closed, piecewise, continuous least-upper-bounding 
constraint surfaces that control the size of successive dependent variable 
excursions arising out of the time-stepping process. 

A list of parameters to be controlled by the self-adaptive solution 
strategies has been generated. The parameters defined to date are listed below 
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Parameters to be Controlled 


1. Element Type(s) 

2. Type(s) of Integration 

3. Order(s) of Numerical Quadrature 

4. Maximum Number of Iterations 

5. Tolerance! s) on Convergence 

6. Constitutive Equation(s) 

7. Yield Criterion (Criteria) 

8. Load Increments 

9. Time Increments 

10. Nonlinear Solution Algorithm(s) 

Experience with the in-house programs has given us a good basis for 
developing the necessary tolerances on convergence. First, convergence is 
evaluated locally, not globally; that is, it is evaluated at each element or 
each numerical integration point. Second, for numerical conditioning, limits 
should be set below which inelastic strains are considered to be zero. Third, 
for time-dependent effects, both temperature and stress cutoffs should be 
estaolished below which time dependent inelastic strain is considered to be 
zero. Then the local convergence criteria for incremental analysis are the 
following . 


Time Independent 


If 


Ep < PCUTOFF , c p 2 0.0 


then 


Aepi < TOL - CONVERGENCE 
or 


4c PI ~ &cp (1-1) 
Acpi 


< TOL - CONVERGENCE 


Time Dependent 


If 


TEMP < TOLC, e c 2 0.0 
and/or 

o e < o c , e c 2 0.0 


87 


then 


A o.j < CTOL - CONVERGENCE 


or 


* q EI Aq E(I-l) 
A°ei 


< CTOL 


CONVERGENCE 


The different convergence criteria are dictated by the wide material strength 
levels encountered in nonlinear analysis. We have also discovered that it is 
advantageous to be able to change these criteria during the course of an 
incremental analysis. 

One approach taken in nonlinear computer codes is the right-hand-side 
technique, in which the plasticity is accounted for by adding an additional 
force vector to the right-hand side of the system of equations 

[K] (d) - (F> ♦ { f p ) 

The basic logic is as follows: 

1. Solve for displacements from 

[K] Cd> - (F } ♦ (f p > 

2. Using the displacements and the constitutive equations, determine ' 
elastic and plastic strains for each element. 

3. Check convergence. 

4. Make an estimate of plastic strains that will satisfy the constitu- 
tive equations, equilibrium, and compatibility. 

5. Based on the estimate of plastic strains from Step 4, form a new 
plastic load vector and go back to Step 1. 

This iteration scheme continues until the convergence criteria are satisfied. 

The plastic iteration accounts for a considerable portion of the total 
computer cost in running a nonlinear finite element code. Substantial improve- 
ments have been made in accelerating the convergence of plastic iteration by 
improving the estimate of the solution in Step 4. Three options are now 
available . 

The first of these schemes is the simplest, and uses the curr*nt calcula- 
tion of plastic strain from the constitutive equations as the estimate of the 
solution. This is the u3ual method on right-hand-side iteration schemes. 
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The second scheme is a modification of the original iteration scheme, and 
is essentially a successive-over-relaxation (SOR) scheme. The estimate of the 
solution is given by: 



t 1 ■ current estimate of solution 
P 


e* - ^ ■ previous calculation of plastic strain from constitutive equations 
P 

e 1 ■ current calculation of plastic strain from constitutive equations 
P 

a = current acceleration factor 

a * 1.5 is used. 

This estimation procedure 
is used: 

c 1 - c 1 ' 1 ♦ 0.5 I' 1 - e 1 * 1 

\ p 

The third scheme is based on an Aitken's extrapolation fonnula for a 
fixed-point iteration. Although we are not reaLly doing a fixed-point itera- 
tion, the finite element equations behave in much the same way. The equation 
used in estimating t{ie solution is: 


:ontinues until e* - < c 1 ^ , then the following 
P P 


~i 

e 


P 


i-2 

e 

P 


i-1 i-2 

e - c 


i 0 T-1 1-2 

e ♦ 2c +e 


Where the symbols are as before and ^ * calculation of plastic strain from 
the constitutive equations two iterations ago. 


The Aitken'-s extrapolation works best when performed every third itera- 
tion. In between Aitken's extrapolations 


a l 

c 

p 


i 

e 


P 


r .oX 
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is used. This equation is also used when the denominator in the Aitken's 
equation approaches zero. 

Test cases using the bolthole model of Figure 72 loaded in tension have 
been run. Comparisons of the number of iterations needed are shown in Figure 
73 for the various iteration schemes. 
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Previous Iteration Scheme 



Figure 73. Effective Plastic Strain in Percent Radius 
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PFRFORM STRUCTURHL nNOLVSlG 





COSMOS DATA DICTIONARY SUMMARY 


DATE 

08/16/83 

08/15/83 

08/16/83 

08/15/83 

08/16/83 

06/20/83 

06/10/83 

06/20/83 

06/07/83 

08/16/83 

07/26/83 

08/16/83 

08/15/83 

08/15/83 

06/10/83 

06/20/83 

08/16/83 

06/10/83 

06/20/83 

06/20/83 

06/13/83 

06/12/83 

08/15/83 

06/19/83 

08/15/83 

06/13/83 

07/26/83 

07/26/83 

07/26/83 

08/16/83 

08/15/83 

08/16/83 


REV NAME 

A2 ANALYSIS CONDITIONS 

AO ANALYSIS INPUT DECK 

A 1 ANALYSIS OUTPUT CONTROL 

AO ANALYSIS RESULTS 

A 2 CONFIGURATION DETAIL GEOMETRY 

A5 CRITICAL POINT HISTORY 

AO CRITICAL POINT LIFE 

A1 CRITICAL POINT LOCATION 

A1 CYCLE MAP 

A2 DECISION PARAMETERS 

AO DECOMPOSED MISSION ELEMENT MATRIX 

A1 DECOMPOSED MISSION ELEMENT MISSION DEFINITION 

A1 DISCRETE MODEL 

AO ELEMENT PRESSURES 

AO ENGINE RATING 

A1 HEAT TRANSFER DATA 

A2 HEAT TRANSFER MATERIAL PROPERTIES 

A 1 HEAT TRANSFER PERFORMANCE DATA 

A1 MATERIAL LIFE DATA 

A 1 METAL TEMPERATURES 

AO MISSION DEFINITION 

A2 MISSION PERFORMANCE DATA 

AO MODEL MODIFICATIONS 

A 2 MODELLING CRITERIA 

AO NODAL TEMPERATURES 

A1 PART GEOMETRY 

AO PART HISTORY 

AO SOLUTION CRITERIA 

AO SOLUTION RESULTS 

A2 STRUCTURAL BOUNDARY CONDITIONS 

A 1 STRUCTURAL DATA 

A2 STRUCTURAL MATERIAL PROPERTIES 


IHI 

»• DATA DICTIONARY ••••• 





DATE: 

08/16/83 



REV: 

A2 

DATA FLOW NAME: ANALYSIS 

CONDITIONS 



ALIASES: NONE 




DIAGRAM REFERENCE: A4 


. 



COMPOSITION: 

Element data: 

Element Number 

♦ {Node Number} 

♦ Element Type 

♦ Material 

♦ Angles for orthotropic properties 

♦ Plasticity indicator 

♦ Creep indicator 

♦ Pressure 
Node Data: 

Node Number 

♦ X f Y,Z location 

♦ Temperature 

♦ Force 
•*- Angles 

♦ Boundary Conditions 


NOTES: 
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DATA DICTIONARY 


DATE: 08/15/83 
REV: AO 


DATA T LOW NAME: ANALYSIS INPUT DECX 

ALIASES: NONE 

DIAGRAM REFERENCE: A4.3 


COMPOSITION: 

ANALYSIS INPUT DECK - 


[Node number] 

+ (Element Number) 

♦ (Temperature) 

♦ (Pressure) 

+ (Boundary Condition) 


NOTES: 


1. The analysis program being used determines the format of 
the Analysis Input Deck. 
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DATA DICTIONARY 


DATE: 08/16/83 
REV: A1 

DATA FLOW NAME: ANALYSIS OUTPUT CONTROL 
ALIASES: NONE 
DIAGRAM REFERENCE: AM 


COMPOSITION: 

ANALYSIS OUTPUT CONTROL s Flag that indicates whether or not 

a Part History should be generated. 


NOTES: 


no 


DATA DICTIONARY 


DATS: 08/15/83 
REV: AO 


DATA FLOW NAME: ANALYSIS RESULTS 

ALIASES; NONE 

DIAGRAM REFERENCE: A4.3 


COMPOSITION: 

ANALYSIS RESULTS 


(Displacement) 

4 (Reaction Force) 
4- (Stress) 

«► (Strain) 


NOTES: 
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DATA DICTIONARY 


DATE: 08/16/83 
REV: A2 

DATA FLOW NAME: CONFIGURATION DETAIL GEOMETRY 
ALIASES: CONFIG DETAIL GEOM 
DIAGRAM REFERENCE: 


COMPOSITION: 


CONFIGURATION DETAIL GEOMETRY = 


({HOLE}) 

(INTERNAL 

(COATING) 


PASSAGE GEOMETRY) 


NOTES: 


1. Part of HEAT TRANSFER DATA 
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•••*• DATA DICTIONARY 


DATE: 06/20/63 
REV: AS 

DATA FLOW NAME: CRITICAL POINT HISTORY 
ALIASES: NONE 
DIAGRAM REFERENCE: AO 


COMPOSITION: 


CRITICAL POINT HISTORY « 


TIME 

♦ STRESS 
< ♦ STRAIN 

♦ DISPLACEMENT 

♦ TEMPERATURE 


NOTES: 


CRITICAL POINT HISTORY contains information about selected 
critical points in a mission. 
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••••• DATA DICTIONARY 


DATE: Ql/'O/ll 
REV: AO 

data~flo£~mame7”critical point life 

ALIASES: NONE 
DIAGRAM REFERENCE: A-0 

COMPOSITION: 

CRITICAL POINT LIFE « Ti*«/cycl«s to "failur#" 


NOTES: 



••••• DATA DICTIONARY ••••• 

DATE: 06/20 '63 

REV: A 1 

DATAFLOW NAME: CRITICAL POINT LOCATION 
ALIASES: CRITICAL POINT LOC 
DIAGRAM REFERENCE: A-0 

COMPOSITION: 

CRITICAL POINT LOCATION * {COORDINATES OF CRITICAL POINT} 


NOTES: 
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DATA DICTIONARY 


DATE: 06/07/83 

REV: A 1 

DATA FLOW NAME: CYCLE MAP 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 

COMPOSITION: 

CYCLE MAP * {CYCLE CASE} 


NOTES: 


1. There are 148 cases in a CYCLE MAP. 
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DATA DICTIONARY 


DATE: 06/16/83 
REV: A2 

DATA FLOW NAME: DECISION PARAMETERS 
ALIASES: DECISION PARAM 
DIAGRAM REFERENCE: A-0 


COMPOSITION: 


DECISION PARAMETERS 


Solution Type 

♦ Structural Analysis Program 

♦ RPM 

♦ Acceleration 

♦ Number of Nodes 

♦ Number of Elements 

♦ Number of different materials 

♦ Number of Pressure Boundary Conditions 

♦ Reference Temperature 

♦ Number of Incremental Load Conditions 

♦ {Constitutive Equation} 

♦ {Output parameter} 

♦ {Re-start parameter} 

♦ {Solution option} 

♦ {Convergence Control} 


NOTES: 

1. Solution Types include: 

Plane Stress 
Plane Strain 
Axisymmetric 
3D 

(with/without Thermal Strain) 
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NOTES : 
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••••• DATA DICTIONARY ••••• 

DATE: 08/16/83 
REV: A1 

DATA FLOW NAME: DECOMPOSED MISSION ELEMENT MISSION DEFINITION 
ALIASES: DME MISSION DEFINITION 
DIAGRAM REFERENCE: ?? 


COMPOSITION: 

DME MISSION DEFINITION 


Altitude 

♦ Mach Number 

♦ Power level 

' ♦ Ambient temperature 

♦ Bleed 

♦ Deterioration 

Time to next point . 


NOTES: 
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DATA DICTIONARY 


DATE: 08/15/83 
REV: A1 


DATA FLOW NAME: DISCRETE MODEL 
ALIASES: NONE 
DIAGRAM REFERENCE: AO 


COMPOSITION: 

PART MODEL - (NODE) 

♦ (ELEMENT) 

♦ (PART GEOMETRY) 


NOTES: 

1. A NODE consists of: 

Node Number 

♦ X coordinate 
+ Y coordinate 
+ Z coordinate 

2. An ELEMENT consists of: 

Element Number 
+ Element Type 

♦ (Node Number) 
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***** DATA DICTIONARY ***** 

DATE* 08/15/83 
REV: AO 

^B ^B mm mm ^b • • « ^m mm ^m ^m ^m ^m « • mm « mm w ^b ^b ^b ^m ^b ^b ^b ^b ^b ^b ^b mm « ^B w » « * ^b «« «• • «b « • «b « «* « w» « ^b ^b ^b ^b ^b « w ^b ^b ^b 

DATA FLOW NAME: ELEMENT PRESSURES 
ALIASES: NONE 
DIAGRAM 'REFERENCE: A4.1 

COMPOSITION: 

ELEMENT PRESSURES - 


NOTES: 
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••••• DATA DICTIONARY 


DATE: 06/10/83 

REV: AO 


DATA FLOW NAME: ENGINE RATING 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 


COMPOSITION: 

ENGINE RATING » Maximum power level at a given altitude, mash 

number, and ambient temperature 


NOTES: 
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•*«•* DATA DICTIONARY •■••• 

DATE: 06/20/83 

REV: A 1 

DATA-FLOW” NAME : HEAT TRANSFER DATA 
ALIASES: HEAT XFER DATA 
DIAGRAM REFERENCE: A-0 

COMPOSITION: 

HEAT TRANSFER DATA * CONFIGURATION DETAIL GEOMETRY 

♦ HEAT TRANSFER MATERIAL PROPERTIES 


NOTES: 
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DATA DICTIONARY ••••• 

DATE: 08/16/83 
REV: A2 

DATA FLOW NAME: HEAT TRANSFER MATERIAL PROPERTIES 
ALIASES: HEAT XFER MATERIAL PROPERTIES 
DIAGRAM REFERENCE: 


COMPOSITION: 

HEAT TRANSFER MATERIAL PROPERTIES * {SPECIFIC HEAT} 

♦ {VISCOSITY} 

♦ {THERMAL CONDUCTIVITY} 

♦ {THERMAL NUMBER} 


NOTES: 


1. Part of HEAT TRANSFER DATA 
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••••• DATA DICTIONARY ••••• 


DATE: 06/10/83 
REV: A 1 


DATA FLOW NAME: HEAT TRANSFER PERFORMANCE DATA 
ALIASES: HEAT XFER PERFORMANCE DATA 
DIAGRAM REFERENCE: AO 


COMPOSITION: 


HEAT TRANSFER PERFORMANCE DATA * 


CASE NUMBER 
RPM 


{ 


GAS PRESSURE 

♦ GAS TEMPERATURE ' 

♦ LOCATION 


NOTES: 

1. Or.c case for each unique mission deflntion point. 
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••••• DATA DICTIONARY 


DATE: 06/19/83 

REV: A 1 

DATAFLOW - NAME: "material LIFE DATA 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 


COMPOSITION : 

MATERIAL DATA « {LOW CYCLE FATIGUE DATA} 

♦ {RUPTURE LIFE DATA) 

♦ {CREEP LIFE DATA) 

♦ {MATERIAL OXIDATION RATE DATA} 


NOTES: 

1. LOW CYCLE FATIGUE DATA: Cyclic life vs stress or strain 

range and metal temperatures. 

2. CREEP/RUPTURE LIFE DATA: Hours to "failure" vs stress and 

■etal temperatures. 

3. MATERIAL OXIDATION RATE DATA: Oxidation rates vs metal 

temperature, gas density, and gas velocity. 
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•*••• DATA DICTIONARY 


DATE: 06/19/83 
REV: A 1 

DATA FLOW NAME: METAL TEMPERATURES 
ALIASES: NONE 
DIAGRAM REFERENCE: AO 


COMPOSITION: 

METAL TEMPERATURES * 


' CASE NUMBER 

, f* TEMPERATURE 1 
+ l ♦ ELEMENT NUMBER } 
[ ♦ NODF. NUMBER J 


NOTES: 
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••••• DATA DICTIONARY 


DATE: 06/13/63 

REV: AO 

DATA FLOW NAME: MISSION DEFINITION 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 


COMPOSITION: 


MISSION DEFINITION s 


' ALTITUDE 

♦ MACH NUMBER 

♦ POWER LEVEL 

- ♦ AMBIENT TEMPERATURE 

♦ BLEED 

♦ DETERIORATION 

.♦ TIME TO NEXT POINT . 


NOTES: 

1. A series of points that define a mission 
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••••• DATA DICTIONARY •«••• 


DATE: 06/12/83 
REV: A2 

DATA FLOW NAME: MISSION PERFORMANCE DATA 
ALIASES: MISSION PERF DATA 
DIAGRAM REFERENCE: AO 


COMPOSITION: 


MISSION PERFORMANCE DATA « 


CASE NUMBER /TIME LIST 


CASE NUMBER 
♦ BPM 



9 


GAS PRESSURE 
♦ LOCATION 


} 




NOTES: 


1. One case for each unique mission definition point. 

2. CASE NUMBER/TIME LIST provides the sequence of CASE NUMBERS 
that define a mission and the time from each case to the 
next. 
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DATA DICTIONARY 


DATE: 08/15/83 
REV: AO 

DATA FLOW NAME: MODEL MODIFICATIONS 
ALIASES: NONE 
DIAGRAM REFERENCE: A4.1 

COMPOSITION: 

MODEL MODIFICATIONS ■ f Node Number / Element Number 

1 + addition-deletion indicator 


NOTES 



DATA DICTIONARY 


DATE: 06/19/83 
REV: A 2 

DATA FLOW NAME: MODELLING CRITERIA 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 

COMPOSITION: 

MODELLING CRITERIA * Parameters used to define discretized mesh 


NOTES: 
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DATA DICTIONARY 


DATS: 08/15/83 
REV : AO 

DATA FLOW NAME: NODAL TEMPERATURES 

ALIASES: NONE 

DIAGRAM REFERENCE: A4 . 1 


COMPOSITION: 

NODAL TEMFERATURES - 


{ 


Node Number I 
+ Temperature | 


NOTES: 


1. Temperature is given in degrees Fahrenheit. 
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DATA DICTIONARY 


DATE: 06/13/83 
REV: A 1 

DATA FLOW NAME: PART GEOMETRY 

ALIASES: NONE 

DIAGRAM REFERENCE: A-0 


COMPOSITION: 

PART GEOMETRY * {POINTS) 

♦ (CURVES) 

♦ (SURFACES) 

♦ (REGIONS) 


NOTES: 
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DATA DICTIONARY 


DATE: 07/26/83 
REV: AO 


DATA FLOW NAME: PART HISTORY 
ALIASES: NONE 
DIAGRAM REFERENCE: A4 

COMPOSITION: 

PART HISTORY * T Stress 
< ♦ Strain 
(_+ Displacement 


NOTES: 
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DATA DICTIONARY 


DATE i 07/26/83 
REV: AO 

DATA FLOW NAME: SOLUTION CRITERIA 
ALIASES: NONE 
DIAGRAM REFERENCE: A4 

COMPOSITION: 

SOLUTION CRITERIA - Soultion Type 

+ 77? 


NOTES 





DATA DICTIONARY 


DATE: 07/26/83 
REV: AO 


DATA PLOW NAME: SOLUTION RESULTS 
ALIASES: NONE 
DIAGRAM REFERENCE: A4 

bomb mb w w a* « « ^B m> « ^B w «m • a • «b«m mb BBaaMasaBwaB mb m w ^b ^b ^b «b ^b «a ^b bb i 

COMPOSITION: 

SOLUTION RESULTS « Maximum stress deviation 

Maximum strain deviation 

NOTES: 
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DATA DICTIONARY 


DATE: 08/16/83 
REV: A2 

DATA FLOW NAME: STRUCTURAL BOUNDARY CONDITIONS 
ALIASES: STRUCT BOUND COND 
DIAGRAM REFERENCE: 


COMPOSITION: 


PART BOUNDARY CONDITIONS 


{DISPLACEMENT} 

{LOAD} 

{STRESS} 

{STIFFNESS} 


NOTES: 


1. . Part of STRUCTURAL DATA 
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DATA DICTIONARY 


DATEs 06/15/83 
REVs A1 


DATA FLOW NAME: STRUCTURAL DATA 

ALIASES s NONE 

DIAGRAM REFERENCES A-0 


COMPOSITIONS 

STRUCTURAL DATA - CONFIGURATION DETAIL GEOMETRY 

♦ (STRUCTURAL BOUNDARY CONDITION} 

♦ (STRUCTURAL MATERIAL PROPERTIES] 
+ {Plasticity indicator) 

+ {Creep indicator} 


NOTES: 
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DATA DICTIONARY 


DATE: 08/16/83 
REV: A2 

DATA FLOW NAME: STRUCTURAL MATERIAL PROPERTIES 
ALIASES: STRUCT MATL PROP 
DIAGRAM REFERENCE: 

COMPOSITION: 

STRUCTURAL MATERIAL PROPERTIES * {STRESS/STRAIN CURVES} 

♦ CREEP DATA 

♦ {Orthotropic angle) 


NOTES: 

1. Part of STRUCTURAL DATA 
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APPENDIX C 

THERMODYNAMIC ENGINE MODEL SPECIFICATION . 
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I. 


II. 


III. 


A. Setup information (furnished with modal) 

1. Engine performance data; 148 cases per Table I-A-l . 

a. Case parameters per Table I-A-l-a. 

2. Engine rating data - Table I -A- 2 . 1-3 . 

3. Power level index matrix, Table I -A- 3 . 

B. User Information 

1. Mission definition data 

a. One line for each mission phase point. 

b. Each line contains the following control data: 

1. Phase # 

2. Mach number 

3. Altitude - feet 

4. Offset from standard day - *F 

5. Power level parameter code # 

6. Power level parameter value 

7. Customer bleed - f/sec. 

8. Deterioration level - *F 

9. Time Increment between this phase point and the 
next, min. 

c. One line, following the mission phase point data line, 
for each parameter to be offset from its steady state 
value. 

d. - Each offset line shall contain: 

1. Phase # 

2. Parameter if 

3. Offset factor 

4. Offset adder 


A. A performance case for each mission phase point 

1. Parameters per Table I-A-l-a. 

2. Format similar Co Table Il-A-2. 


Technical Basis 

A. Each new case will be generated from available cases (I-A-l) by 
a disciplined interpolation process similar to that currently 
used in the Life Analysis by Stress and Temperature Simulation 
(LASTS) program. 

1. All parameters will be transformed to a functional form that 
has optimal linearity relative to all other parameters, 
a. A study will be performed on CF6-50C2 engine 
performance data to evaluate and improve the 
interpolation functions. 
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2. Each transformed parameter vlll b« interpolated by a 
linear interpolation procaaa and chan cransforaad back 
co lea noraal fora. 

3. Tha lncarpolacion "targets" ahall ba apaclflad in Che 
lnpuc for aach mission phaaa point (I-B-l-b). 

4. Tha lncarpolacion procss* for aach phaaa point ahall begin 
with a base caaa near .he desired mission phase point. 

5. Each lncarpolacion step will convert cha base case (or 
previously modified case) to tha desired level of Cha 
target control parameter (l.e., MN) . 

a. Linear partlals are assumed; Interactions are ignored. 

b. Partlals vlll be derived from tvo or more "partial 
cases" near ehe base case conditions. 

c. The interpolation steps will be performed in Che 
sequence in Table I-A-l-s. 

d. The specific power level parameter to.be used as an 
interpolation carget la input as "power level #", 
followed by the target value. 

e. For flight idle and ground idle, a special power 
level # (one for FI or one for Cl) will be entered, 
followed by a zero parameter value; the standard 

FI and GX power levels will be used. 

f. For thrust reverse, a special power level # will be 
used and a value of fan speed vlll follow; thus 
thrust reverse power level will always be based on a 
fan speed target. 

The LASTS interpolation process will be modified to eliminate 
the current manual procedure for the generation of interpolation 
instructions. 

1. A sec of base case numbers will be provided as a function of 
altitude and Mach number. * 

2. A family of pairs or triplets of partials case numbers will 
be provided as a function of altitude and Mach number for 
each control parameter. 

3. The user will be required to input only the data in I-B-l-b 
for each mission phase point. 

The Thermodynamic Engine Model (TDE) shall have the capability 
to predict the minimum time for speed changes due to throttle 
actions. 

1. The user shall have the option to input zero throttle-action 
transient times, and the model will calculate appropriate 
transient times, subtracting them from following phase 
times . 

2. The transient time calculation shall be sensitive co the 
effects of altitude. 

The user shall have the option of selecting an appropriate 
CF6-50C2 power management point, and avoid the need for 
specifying absolute values of the power level parameter. 

1. The power management parameter code shall call for take-eff, 
max climb, or max cruise rating, and the adjacent value 
shall specify the I derate desired based on Z thrust. 


E. Offsets of specific parameters, relative to the steady state 
performance cases shall be permitted to simulate take-off 
transient conditions. 

1. Each parameter change shall be specified by a line 
following che mission phase point data. 

2. The parameter offset data shall be: 


a. 

b. 

Case # 
Parameter # 

c. 

Offset 

factor 

d. 

Offset 

adder 


3. Offset calculations vlll be performed after mission phase 
point Interpolations are completed. 

Software Character is tics /Interfaces 

A. Later (to be integrated with overall software of the COSMOS 
Program) . 
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CF6-50C2 PERFORMANCE CASES 
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TABLE I-A-l 

CF6-50C2 PERFORMANCE CASES 
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TABLE I-A-l 

CF6-50C2 PERFORMANCE CASES 


M 
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TABLE I -A- la 

ENGINE PERFORMANCE CASE PARAMETERS 


P2 

Fan Inlet Total Pressure 

PSIA 

P3 

Compressor Discharge Total Pressure 

PSIA 

P4 

Turbine Inlet Total Pressure 

PSIA 

P49 

Turbine Outlet Total Pressure 

PSIA 

T2 

Fan Inlet Total Temperature 

•f 

T3 

Compressor Discharge Total Temperature 

*F 

T41 

Turbine Inlet Total Temperature 

•f 

T49 

Turbine Outlet Total Temperature 

•f 

W25 

Fan Air Flow 

#/ sec 

FNIN1 

Installed Thrust 

t 

DTAMB 

Offset from Standard Day Temperature 

•f 

W41 

Turbine Air Flow 

#/sec 

XNB 

Core Speed 

RPM 

XNL 

Fan Speed 

RPM 

MN 

Mach Number 


ALT 

Altitude 

Feet 

WB27/VB3 

Customer Bleed 

#/ sec 

DT49 

Engine Deterioration Index 

°F 
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TABLE I-A-2.1 

APPROX. CF6-50C2 RATING DATA 
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TABLE I-A-2.2 

APPROX. CF6-50C2 RATING DATA 
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21.2 103.4 3549 

34.0 98.5 3381 98.5 3381 

60.0 90.4 3103 90.4 3103 90.4 3103 


TABLE I-A-2.3 

APPROX. CF6-50C2 RATING DATA 
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20.0 

20.2 

<\SJ 

bO.O 07. 35 


TABLE I-A-3 


AUTOMATED INTERPOLATION SYSTEM 
POWER LEVEL INDEX 
CF6-50C2 

' PCNLR 



SET 

SET 

P.L. INDEX 

0 

10000’ 

1 

109.0 

119.0 

2 

100.0 

109.0 

3 

90.0 

98.0 

4 

75.0 

81 .0 

5 

55.0 

53.0 

6 

38.0 

44.0 

7 

24.0 

34.0 


ALT 

ALT 

P.L. INDEX 

10000’ 

40000’ 

1 

103.5 

120.0 

2 

91.5 

106.0 

3 

75.0 

85.0 

4 

58.0 

67.0 

5 

44.0 

48.0 

6 

34.0 

37.0 
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original page ,3 
0F POOR quality 


TABLE II-A-2 

PERFORMANCE OATA OUTPUT FORMAT 


• ENGINE PERFORMANCE DATA BY MISSION PHASE ■ 


phase 

» 

i 

CASE 

0 

1*7 

•• E N G I 
P48 

14.476 

N E PAR 

P3 

14.479 

A M E T E R •• 

V2S FNIN1 OTAMB 

1.108 0. 1.700 

2 

190 

14.476 

14.478 

1. 108 

0 . 

1.700 

3 

1*1 

14.47* 

14.476 

1.106 

0 . 

1.700 

4 

192 

13. 179 

14.476 

1. 108 

0 . 

1.700 

9 

193 

16.476 

16.476 

1. 108 

0 . 

1.700 

« 

131 

19.912 

14.393 

38 . 302 

1943.746 

1.700 

7 

131 

16.912 

14.393 

38.302 

1943.746 

1.700 

• 

132 

80.731 

20 . 930 

243.413 

42726.324 

1.701 

9 

133 

82.833 

21.307 

243.733 

39699.800 

2.301 

10 

134 

66. 172 

18.773 

206 . 290 

30340.937 

2.302 

1 1 

133 

63.487 

16.303 

194.736 

21431.468 

3.001 

12 

136 

34.321 

8. 196 

103. 120 

9766.340 

12.200 

13 

137 

24 . 340 

3.313 

81. 138 

6421.238 

12.200 

14 

137 

24.340 

3.315 

81.138 

6421.258 

12.200 

13 

138 

4.388 

2.7*7 

19.617 

-892.969 

12.200 

1« 

139 

8.911 

6.888 

26.739 

- 1914.288 

7.800 

17 

160 

13.099 

10.748 

33 . 303 

-2310.870 

4.200 

18 

161 

16.297 

13.944 

38. 100 

894. ns 

2.500 

19 

162 

27.946 

14.707 

93 . 330 

7273.440 

2.500 

20 

163 

29. 131 

13.402 

99.313 

8873.031 

1.700 

21 

164 

22.015 

14.861 

68.412 

4493. 193 

1.700 

22 

164 

22.013 

14.861 

68.4 i2 

4493. >93 

1.700 

23 

163 

43.788 

16.292 

143.437 

14*71. 125 

1.7C4 

24 

166 

43.084 

16. 184 

144.799 

13603.446 

1 700 

23 

187 

16.942 

14.822 

38.538 

1946.327 

1.700 

2« 

187 

16.942 

14.622 

38.358 

1946.327 

1.700 

27 

194 

16.303 

16.303 

0.974 

0 . 

1.700 

28 

193 

14.503 

14.303 

0.974 

0 . 

1.700 

29 

196 

14.303 

14.303 

0.974 

0 . 

1 . 700 

30 

186 

14.303 

14.303 

0.974 

0 . 

1.700 
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